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We introduce a concise formalism to characterize nanometer-sized tori based on carbon nanotubes 
and to determine their stability by combining ah initio density functional calculations with a contin¬ 
uum elasticity theory approach that requires only shape information. We find that the high strain 
energy in nanotori containing only hexagonal rings is significantly reduced in nanotori containing 
also other polygons. Our approach allows to determine local curvature and link it to local strain 
energy, which is correlated with local stability and chemical reactivity. 
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INTRODUCTION 

The past three decades witnessed an unprecedented 
progress in the synthesis and use of nanostructures con¬ 
sisting of elemental carbon ranging from fullerenes to 
carbon nanotubes and graphene^^^. Due to the stiff¬ 
ness of the (j bonding network and the delocalization of 
the TT electrons in these systems, the structure-property 
relationships are highly nontrivial, as presence of local 
atomic defects is known to affect physical properties on 
the mesoscopic length scale^. In this regard, owing to 
their peculiar geometry and topology^, carbon nanotori 
based on single-wall carbon nanotubes serve as an excel¬ 
lent test ground for studying such effects. To name just 
a few, reversible elastic deformation of carbon nanotori 
has been predicted^ and experimentally identified^, im¬ 
plying potential application in force-sensing devices. Gi¬ 
gantic paramagnetic response^, persistent currents^^”^^ 
and even molecular anapole moments^^ have been pre¬ 
dicted for specific nanotorus geometries. Also, the in¬ 
tricate shape of carbon nanotori can be exploited in the 
context of host-guest chemistry and physics; for exam¬ 
ple, megahertz oscillations of a nanotorus mounted on 
a nanotube^^ and metal-encapsulated nanotori with net 
magnetic momenthave been predicted. A crucial pre¬ 
requisite to realize these interesting physical properties 
is to establish, to what degree the postulated toroidal 
nanostructures are stable mechanically and thermally. 

Here we study the stability of carbon nanotori formed 
conceptually by bending a finite-length single-wall car¬ 
bon nanotube (GNT) and connecting its ends seamlessly. 
Such a nanotorus can then be characterized the chiral 
index^^ (n, m) of the nanotube and number of primitive 
unit cells along the perimeter. We will refer to this fam¬ 
ily of nanotori, which contain only hexagonal rings, as 
to polyhex nanotori. It is known that presence of non- 
hexagonal rings, such as 5-7 pairs, induces a bend in a 
straight nanotube. In the following, we will refer to the 
family of nanotori that contain also non-hexagonal rings 
as to polygonal nanotori. We will demonstrate that a con¬ 
tinuum elasticity theory approach^^, which requires only 


shape information, is capable of estimating not only the 
global stability, but also local strain with a precision ap¬ 
proaching that of ab initio density functional calculations 
at a small fraction of the computational cost. 

Our manuscript is organized as follows. In the next sec¬ 
tion we outline the computational techniques used in our 
study. Next we summarize the characteristics of different 
nanotori families. Then, we analyze and obtain analyti¬ 
cal expressions for the in-plane and out-of-plane contri¬ 
butions to the elastic energy of polyhex nanotori. In the 
following Section, we will present a general overview of 
the elastic energy for polygonal nanotori. Next we will 
discuss the distributions of curvature and elastic energy 
on polygonal nanotori, and correlate them to the relative 
positions of the non-hexagonal rings in the systems. This 
is followed by the asymptotic analysis of three families of 
polygonal nanotori with varying shape parameters, in¬ 
cluding changing rotational symmetry, the lateral torus 
size normal to the axis and the torus height along the 
axis. 


COMPUTATIONAL TECHNIQUES 

Our main computational approach is based on contin¬ 
uum elasticity theory, which proved useful in evaluating 
the strain energy in fullerenes, nanotubes and schwarzites 
with respect to a planar graphene monolayer and corre¬ 
sponding stability differences^^. 

In systems with no frustration, where areas of Voronoi 
polygons associated with individual atoms are at the op¬ 
timum value, the in-plane strain energy is very small. 
This proved to be the case in carbon nanotubes^^ and 
we expect it to hold also in the related nanotori. Ne¬ 
glecting the in-plane component of strain, the stability of 
nanotori could be analyzed in the inextensional limit of 
the mechanical theory of a membrane^^, where the elastic 
energy is dominated by its out-of-plane component. 

The out-of-plane strain or curvature energy of a given 



2 


closed membrane can be expressed as 

AEc = D [ dA[2k‘^ - (1 - a)G] , (1) 

Js 

where D is the flexural rigidity, a the Poisson ratio, k 
the local mean curvature, G the local Gaussian curva¬ 
ture, and the integration is carried out over the surface 
S of the membrane. If instead of the general shape we 
know the precise atomic locations from either diffraction 
experiments or ab initio calculations, we can rather use 
the discretized version of Eq. (1), 

AEc « AqD - (1 - a)Gi] . (2) 

i 

Here, Aq is the area per atom, ki the local mean curva¬ 
ture and Gi the local Gaussian curvature at atom site i. 
The summation covers all atoms of the torus. For struc¬ 
tures comprising of more than one element, the quantity 
Aq may vary and needs to be taken inside of the sum¬ 
mation. Alongside with this expression, accurate ways 
to estimate ki and Gi based on the local geometry in the 
vicinity of site i have been proposed recentlywhich 
have been validated by ab initio density functional cal¬ 
culations. Even though this method relies on the knowl¬ 
edge of the molecular geometry, it has been shown that 
even inexpensive classical force fields serve the purpose of 
structure optimization surprisingly well when compared 
against state-of-the-art ab initio calculations for a wide 
range of graphitic structures'^’ 

The total in-plane strain energy can be estimated in 
the continuum model from an integral over the entire 
surface area A, 

AEs = ^ e{a)dA , (3) 

where e(cr) is the energy cost per area subject to strain 
a. For specific, small strains, e(cr) may be presented as 
a harmonic function of the strain cr, with the proper co¬ 
efficient taken either from experimental data or form ab 
initio calculations. 

Selected results based on continuum elasticity the¬ 
ory are validated by ab initio calculations in the frame¬ 
work of density functional theory (DFT), as imple¬ 
mented in the SIESTA code.^^ We use the local den¬ 
sity approximation^to describe exchange and corre¬ 
lation in the system, norm-conserving Troullier-Martins 
pseudopotentials^^, and a double-^ basis including polar¬ 
ization orbitals. In the reference system, we sample the 
2D Brillouin zone of graphene by 16x16 /c-points.^^ The 
small Brillouin zones of isolated nanotori are sampled by 
only 1 /c-point. We use a mesh cutoff energy of 180 Ry to 
determine the self-consistent charge density, which pro¬ 
vides precision in total energy of <2 meV/atom. All 
geometries are optimized using the conjugate gradient 
method,until none of the residual Hellmann-Feynman 
forces exceeds 10“^ eV/A. 



FIG. 1. (Color online) Atomic structure of (a) a polyhex 
nanotorus with only hexagons and (b) a polygonal nanotorus 
containing non-hexagonal rings, (c) A model of a perfect torus 
with major radius R and minor radius r. 0 is the longitudinal 
angle and cj) is the azimuthal angle. The yellow shading in (a) 
and (b) indicates specific polygons in the nanotori. 


NANOTORUS CHARACTERISTICS 

Due to the extreme flexural stiffness of straight 
CNTs 26’27, it is expected that the number of carbon 
atoms of a stable polyhex nanotorus should be at least 
on the order of 10^. An example of a polyhex nanotorus 
containing 400 carbon atoms is shown in Fig. 1(a). For 
the two shaded hexagons located in the inner and the 
outer rims of the nanotorus, it is obvious that the in¬ 
plane strain energy is inversely proportional to the size 
of the nanotori. This will be addressed in detail in the 
following section. 

In contrast to polyhex nanotori, there has been consid¬ 
erable effort in describing the second family of polygonal 
nanotori, which also contain non-hexagonal rings. Due 
to the requirement of trivalency of sp^ carbon atoms and 
the Euler’s theorem, the number of non-hexagonal rings 
on a general closed surface is constrained by 

^ (6 - m)Nm = 6x, (4) 

m>2 

where Nm is the number of the m-gonal rings and y is 
the Euler characteristic of the closed surface in question. 
For a torus, we have y = 0. Thus, if we restrict ourselves 
to polygonal nanotori having only pentagons, hexagons, 
and heptagons, the number of pentagons (A^s) must equal 
that of heptagons (Nj). Apart from this simple equa¬ 
tion, there are many more combinatorial and geometri¬ 
cal constraints for a stable polygonal carbon nanotorus 
to exist^^. The polygonal nanotori covered in this study 
are restricted to those of high symmetry (D^^ or D^/j, 
point group) and a smaller number of atoms (< 60 atoms 
per rotational unit cell). One such polygonal nanotorus 
is shown in Fig. 1(b), a representative pentagon on the 
outer and a heptagon on the inner side of the nanotorus 
are highlighted. As suggested before, the existence of 
these non-hexagonal rings induces natural bending and 
thus drastically reduces the in-plane strain energy. 


ELASTIC ENERGY IN NANOTORI 

As suggested above, the elastic energy within carbon 
nanotori can be divided into two parts: in-plane and 
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FIG. 2. (Color online) (a) Schematic of a striped nanotube 
(left), which is bent to a nanotorus (right), (b) Contour plot 
representing the in-plane strain energy AEs{R,r) according 
to Eq. (9). 


out-of-plane strain energy. The in-plane strain originates 
from an uneven distribution of atomic area densities as 
compared to graphene at equilibrium. This is most obvi¬ 
ous in polyhex nanotori, where hexagons along the outer 
perimeter are stretched and those along the inner perime¬ 
ter are compressed, as seen in Fig. 1(a). The out-of-plane 
strain energy is associated with forming any curved sur¬ 
face of an intrinsically flat material. 

In this section we examine these two contributions to 
the elastic energy of nanotori. In particular, we will pro¬ 
vide and analyze analytical expressions for polyhex nan¬ 
otori without buckling^^. For polygonal nanotori, we will 
show that the contribution from the in-plane strain is 
negligible, similar to the case of small fullerenes^^. 


In-Plane Strain Energy in Polyhex Nanotori 

Points on the surface of a perfect parametric torus are 
described by 

X = {R — r cos 0) cos 0 , 

< y = {R — r cos 0) sin (j) , (5) 

z = r sin 0 . 

\ 

Here and in Fig. 1(c), R is the major and r the mi¬ 
nor radius, 0 the zenith angle, and (j) the azimuthal an¬ 
gle. Note that R > r is required for a normal ring torus 
without self-intersection. A polyhex nanotorus, shown 
in Fig. 1(a), can be constructed by rolling up a finite 
CNT and connecting the two ends together, as indicated 
in Fig. 2(a). The major radius R of the polyhex nan¬ 
otorus can be approximated by the length of the original 
nanotube divided by 27r. In this case, the area density of 
atoms is smaller than the average along the outer perime¬ 
ter of the torus, characterized by 7r/2 < 0 < 37r/2, and 


larger along the inner perimeter. Obviously, the in-plane 
strain will be negligible in comparison to the out-of-plane 
strain if >> r. Should this not be the case, then the in¬ 
plane strain energy will be important, as it describes the 
significant distortion of hexagonal rings along the inner 
and the outer perimeters of a poly hex nanotorus, which 
are highlighted in Fig. 1(a). 

Here we present a simple estimation of the in-plane 
strain energy for polyhex nanotori without buckling. As 
shown in Fig. 2(a), a CNT can be seen as parallel nar¬ 
row strips or nanoribbons of equal length 27vR that are 
connected side-by-side. After the CNT is deformed to a 
nanotorus, the strips on the inside of the nanotorus are 
compressed and the strips on the outside are stretched. 
We further assume that the contributions from all other 
in-plane distortion modes are negligible. As shown in 
Fig. 2(a), the width of the strips is rdO and the length is 
27r{R — r cos 0). The strain a of the strips can be written 
as 


(7(6') = 


27r{R — r cos 0) — 2t:R 

2^r^ 


r 

— — cos 0 . 
ft 


(6) 


The total in-plane strain energy of the nanotorus is just 
the sum of the strain energies of all the strips, 


AE, = 


1 ' 


e(cr)27r(i? — r cos 0)rd0 . 


(7) 


Here, e(cr) is the energy cost per area of graphene subject 
to uniaxial strain a. For sufficiently small strain (cr<5%), 
e(cr) is a parabolic function of a, 

e(cr) = ca(0)‘^ = c-^ cos^ 0 . (8) 

R^ 

A fit to graphene data yields the numerical value c = 
9.94 eV/A^. Substituting Eq. (8) into Eq. (7), we obtain 
an expression for the total in-plane strain energy of a 
nanotorus 

^27r p 

AEs = / cos^ 02ti{R — r cos 0)rd0 

Jo 

= 277 ^ 4 - ( 9 ) 

The behavior of AEs{R^r) in polyhex nanotori is shown 
in Fig. 2(b). As suggested by Eq. (9), the strain energy is 
proportional to and inversely proportional to R. Note 
that this expression is obtained in the harmonic limit of 
small strain, which translates into R » r. Also, we 
have assumed that the nanotori are ideal and described 
by Eq. (5). In reality, the cross-section of an elastic nan¬ 
otorus may deviate from a perfect circle upon relaxation, 
and buckling may occur along the perimeter. In spite of 
its limits, this expression provides a concrete estimate of 
the strain energy of experimentally observed rings based 
on CNTs^^“^^. In the reported cases, R is of the order of 
microns and r of the order of nanometers, so all assump¬ 
tions in deriving Eq. (9) are justified. 
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FIG. 3. (Color online) Definition of the principal radius R 2 of a nanotorus in (a) its interior and (b) its exterior part as a 
function of the zenith angle 0. (c) Contour plot displaying the out-of-plane curvature energy AEc{R,r) of nanotori according 
to Eq. (11). 


Out-of-Plane Strain Energy in a Perfect Parametric 
Torus 


Besides the in-plane strain, there is also an out-of-plane 
strain caused by deviation from planarity, which is rep¬ 
resented by Eq. (1). For the parametric torus given in 
Eq. (5), we infer the infinitesimal area dA in Eq. (1) to be 
dA = rdO '{R — r cos ^)d0. Obviously r is one of the prin¬ 
cipal radii of curvature, which is arbitrarily chosen to be 
Ri. The situation is more complicated for i? 2 , the other 
principal radius of curvature. As illustrated in Fig. 3, R 2 
is given by 


R2{0) = 


R sec 0 — r for — 


r — RsecO 


I <^<f 

Stt 


for + I < 6» < ^ . 


( 10 ) 


Clearly, R 2 ^ oc when 6> = ± J, which is indeed the case 
since the torus is tangent to the planes z = ±r. 

Now the evaluation of AEc according to Eq. (1) is 
straightforward. Since the integrand is independent of 
0, we end up with the expression 

p27T 

AEc{R,r) = irDr / dO{R —r cos 0)x 

Jo 


1 


1 


+ 


2a 


= 27r^D 


Rj R2W " RMO) 

R^ 


ry^{R^r){R-r) 


( 11 ) 


In contrast to perfectly spherical fullerenes, where AEc 
is independent of radius^^, the curvature energy for nan¬ 
otori depends both on the major radius R and the minor 
radius r. The out-of-plane curvature energy AEc{R^r) 
of carbon nanotori is presented in Fig. 3(c). 

We find that the curvature energy does not depend 
on the Poisson ratio a, a result from the Gauss-Bonnet 
theorem which states the Gaussian curvature integrated 
over a closed surface equals the Euler characteristic (x) 
of the surface times 27r. Since x = 0 for nanotori or any 
tubular structures, AEc does not depend on a. Note 
that the IJVR — r term in the expression diverges at 
R = r, corresponding to a horn torus, and becomes imag¬ 
inary for < r, representing spindle nanotori with self¬ 


intersection in the center. In addition, taking the limit of 
^ 00 , this expression reduces to the curvature energy 
of a straight cylinder of length I = 27ri?, 


AEc 


7tDI 

r 


( 12 ) 


Also, by setting the derivative of AEc with respect to 
R to zero, one obtains the minimal torus curvature en¬ 
ergy of AEc = 47r^T) for R = v^r, as also reported 
previously^^’^^. 

As mentioned before, for polyhex nanotori described 
by Eq. (5) and R » r, the in-plane strain energy can 
be adequately quantified by Eq. (9) and the out-of-plane 
strain energy by Eq. (12). Taking polyhex nanotori con¬ 
structed from (10,10) CNTs with r = 7 A as an exam¬ 
ple, the constraint on the major radius R for the har¬ 
monic approximation in Eq. (8) to hold is estimated to 
be i?>150 A. For R = 150 A we find that the in-plane 
strain energy AEg = 450 eV and the out-of-plane strain 
energy AEc = 600 eV are of the same order. On the other 
hand, the majority of rings^^ synthesized by ultrasonic 
wave treatment of CNTs in solution has R = 3500 A. 
In those rings, the contribution from the in-plane strain, 
AEg = 20 eV, is negligible in comparison with the out- 
of-plane strain of AEc = 1.4x10^ eV. 


Strain Energy in Polygonal Nanotori 

In general, the shape of a polygonal nanotorus deviates 
from that of a perfect parametric torus due to the pres¬ 
ence of non-hexagonal defects in a hexagonal network. 
The largest changes in the curvature occur near individ¬ 
ual non-hexagonal defects, and often there are relatively 
flat segments in-between the defects. Consequently, the 
analytical expression in Eq. (11) should only be used as 
a very approximate way to estimate the strain energy of 
a polygonal nanotorus that is roughly characterized by a 
set of effective torus radii (i?, r). For a quantitatively bet¬ 
ter energy estimate, we must use the discretized version 
of Eq. (1), Eq. (2), that takes into account the specific 
shape of a carbon nanotorus. As mentioned previously. 
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we first obtain the molecular geometry from the classical 
Keating force field^^. 

Once the optimum geometry is established, the local 
mean curvature ki and the local Gaussian curvature Gi 
are determined everywhere and substituted into Eq. (2). 
The results of this methodology to a subset of 22 polyg¬ 
onal nanotori are shown in Fig. 4(a). Here the strain en¬ 
ergy calculated through Eq. (2) with geometry optimized 
by Keating potential is represented by the green squares, 
energy calculated with the accurate DFT method by 
black dots, and the Keating potential energy for the 
Keating-optimized optimized geometry by red rhombi. 
Our results show clearly that for Keating-optimized ge¬ 
ometries, strain energies based on Eq. (2) reproduce our 
ah initio results rather well. On the other side, strain en¬ 
ergies estimated using Keating potential alone not only 
significantly underestimates the strain, but also do not 
follow the correct general trend. This firmly establishes 
the applicability of the continuum methodology to polyg¬ 
onal nanotori under investigation. 

For specific nanotori, local geometric features that are 
not described by Eq. (5) contribute to strain energy in 
addition to Eq. (11). Corresponding results are presented 
in Fig. 4(b). We considered a large set of polygonal nan¬ 
otori, which were optimized by the Keating potential, 
and fitted a pair of torus radii (i?, r) for each of them. 
We then correlated the strain energy obtained using 
the analytical expression in Eq. (11) with the more proper 
value based on the optimum discrete geome¬ 

try, based on Eq. (2). As expected, most data points lie 
below the dashed line. Even though 

the strain energies estimated using the two approaches 
appear proportional to each other, it is clear that the an¬ 
alytical expression in Eq. (11) underestimates strain in 
polygonal nanotori significantly. We conclude that while 
a quick estimate of the strain energy based on Eq. (11) 





Local curvature energy AE^/A 



FIG. 5. (Color online)Local Gaussian curvature G (left pan¬ 
els) and local curvature energy AEcjA (right panels) across 
the surface of selected nanotori. Representative examples 
shown are (a) the flattened C 320 nanotorus with 320 atoms 
and (b) the elongated C 280 nanotorus with 280 atoms. G and 
AEcjA are interpolated from their values at the atomic sites. 


is useful once the global parameters (i?, r) are known, 
Eq. (2) should be used if quantitative comparison among 
different structures is intended. 


LOCAL STABILITY ANALYSIS OF POLYGONAL 
NANOTORI 

In the previous Section, we have considered nanotori 
as whole objects and established elasticity theory as a 
valid tool to estimate strain energy with respect to pla¬ 
nar graphene in different torus isomers. The total strain 
energy AE, according to Eqs. (1),(2), (7), may be repre¬ 
sented as an integral or a sum of local contributions. In 
other words, unlike in more complex DFT calculations, 
this approach allows to estimate local contributions to¬ 
wards the total strain energy. In the following, we will 
investigate the local strain distributions within individual 
nanotori. Next, we will apply this approach to analyze 
three series of polygonal nanotori and discuss trends in 
the strain energy in the light of a corresponding asymp¬ 
totic analysis. 


Local Curvature Distribution 


FIG. 4. (Color online) Strain energy AE in different polygo¬ 
nal nanotori, (a) Energy differences per atom AE^f'^ based 
on DFT and based on the Keating potential are 

presented next to curvature energies based on 

Eq. (2) for Keating-optimized geometries. See the main text 
for detailed description, (b) Comparison between the ana¬ 
lytical expression in Eq. (11) for the curvature energy AE^ 
based on the elastic description of a parametric torus and 
the curvature energy obtained using Eq. (2) for 

Keating-optimized discrete torus geometries. 


The local geometry of a smooth, two-dimensional ob¬ 
ject may be characterized by two independent quantities, 
such as the local mean curvature k and the local Gaus¬ 
sian curvature G. These quantities can be used in Eq. (1) 
or its discretized counterpart, Eq. (2), to determine the 
curvature energy. The distribution of the local Gaussian 
curvature Gi and the local curvature energy per area, 
AEc^/A q = D[2k‘l — (1 — (y)Gi], across the surface of 
two representative nanotori is shown in Fig. 5. 

In both cases, and others shown in the Supplemental 
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Material^^, the positively curved segments, shown in red 
in the left panels, are concentrated near the loci of pen¬ 
tagons along the outer perimeter. The negatively curved 
segments, shown in blue, are concentrated near the loci 
of heptagons along the inner perimeter. Specifically, the 
heptagons along the inner perimeter of the C 320 nan¬ 
otorus in Fig. 5(a) are well separated from the pentagons 
along the outer perimeter. This is different from the ax¬ 
ially elongated C 280 torus of Fig. 5(b), where pentagon- 
heptagon pairs a the upper and lower planes form an 
azulene-like pattern. As a consequence, the Gaussian 
curvature is more evenly distributed in the latter. We 
emphasize again that for closed nanotori, the summation 
of the local Gaussian curvatures is strictly zero as dic¬ 
tated by the Gauss-Bonnet theorem. 


The distribution of the local curvature energy is even 
more intriguing. Within many nanotori we investi¬ 
gated, some of which are discussed in the Supplemental 
Material^^, we observed some degree of correlation be¬ 
tween the absolute value of the local Gaussian curvature 
and the curvature energy. This should imply, at least for 
the nanotori investigated, that the two local curvatures 
ki and Gi are not entirely independent. Yet the worth 
of this correlation has its limitations, as shown in Fig. 5. 
Whereas in the flattened G320 nanotorus in Fig. 5(a) the 
curvature energy is rather evenly distributed across the 
structure, the strain is clearly largest near the upper and 
lower ends of the G 280 nanotorus in Fig. 5(b). This cur¬ 
vature energy distribution in the right panels differs ob¬ 
viously from the Gaussian curvature distribution in the 
left panels. The reason for this finding is that in these 
extreme structures, we can not truly decouple ki and Gi. 
In G 28 O 5 and ki is constant in the central Tubular 

segments’. Only at the upper and lower ends, a large 
mean curvature ki is required to connect the inner and 
the outer tube. The flatter G320 nanotorus lacks Tubu¬ 
lar segments’ with Gi^O. Therefore, the Gaussian curva¬ 
ture and curvature energy are better correlated and more 
evenly distributed in this isomer. The different shapes of 
carbon nanotori will be discussed in more detail later on. 


We need to point out that the determination of the 
local Gaussian curvature Gi and the local mean curva¬ 
ture ki is a non-trivial task that requires extra attention 
in discrete, irregular structures. According to the proce¬ 
dure outlined in Ref. 17, the determination of Gi requires 
second-nearest-neighbor information. On the other hand, 
the trivalency of graphitic carbon system lends itself to a 
compact definition of ki based on nearest neighbors only. 
The asymmetry in handling the two curvatures leads to 
the possibility of negative local curvature energies accord¬ 
ing to Eq. (2), which is unphysical. This can be resolved 
by lowering the resolution of ki^ e.g. averaging its values 
using a window that also includes second-nearest neigh¬ 
bors. 



FIG. 6. (Color online) (a) A schematic model of a polygo¬ 
nal nanotorus with n-fold symmetry (n = 6 is shown here), 
consisting of n nanotube segments that are connected by n 
nanotube elbow joints. The nanotube segments are character¬ 
ized by the width VF, length L and height H. (b-d) Curvature 
energies per atom based on Eq. (2) for Keating-optimized nan¬ 
otori of different shapes. AEc is presented for nanotori with 
(b) n,W,H = const., but nanotube length L changing, (c) 
n,L,W = const., but nanotube height H changing, and (d) 
L,W,H = const., but rotational symmetry number n chang¬ 
ing. The insets show nanotorus models with different values 
of L, if and n. The lines in (b) and (c) are analytical ex¬ 
trapolations discussed in the text. Two sets of data points in 
(d), connected by lines to guide the eye, correspond to two 
structural families of nanotori, described in the text. 


Shape Dependence of Stability 

In spite of a large number of published 
reports^^’^^’^^”^^, a systematic investigation of how 
the nanotorus stability depends on its shape has been 
missing. Here we intend to cover this gap by examining 
the dependence of the total curvature energy on main 
parameters defining the shape of polygonal carbon 
nanotori, namely the rotational symmetry number 
n, the length L and the width W of the nanotorus 
segments, and the height H of the nanotorus, shown 
schematically in Eig. 6(a). A more detailed classification 
scheme of polygonal nanotori is provided in References 
28 and 47. In the following, we discuss the dependence 
of the curvature energy per atom AEc on the shape 
parameters. The structures used to obtain the numerical 
results in Eig. 6(b-d) are displayed in the Supplemental 
Material^^, along with the distribution of local Gaussian 
curvature and curvature energy following the convention 
used in Eig. 5. 


Side Length Dependenee 

We first look at an interesting case of a polygonal nan¬ 
otorus consisting of n straight nanotube segments of fi¬ 
nite length L that are connected by n elbow joints as 
illustrated in Eig. 6(a). As suggested in the insets of 
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Fig. 6(b), increasing the length of the nanotube segments 
reduces gradually the influence of the joints and for large 
L, the energy of the elbow joints becomes negligibly small 
when compared to the curvature energy of the straight 
tubular segments. This is illustrated for a specific family 
of nanotori in Fig. 6(b), where the data points repre¬ 
sent results obtained using Eq. (2) for a set of Keating- 
optimized nanotori with different side length L that are 
displayed in the Supplemental Material^^. 

For this particular series of nanotori, we may express 
the curvature energy per atom by AEc{L)/N = a -f b/L 
and display this dependence, with a = 0.254 eV and b = 
0.357 eV-A, by the dashed line in Fig. 6(b). 

A more specific expression can be derived assuming 
that the nanotube segments are characterized by the 
chiral index (m,n). In that case, the nanotube radius 
r = 1.42 A X^/3/(27r) can be used in 
Eq. (12). Eor a nanotube segment of length L, the num¬ 
ber of atoms can be estimated using N = 2'KrLjA{^^ 
where Aq = (1.42 A)^x3v^/4 = 2.62 A^ is the area 
per atom in graphene. Then, we obtain 

AEc{m,n) ^ ^ AEc ^ V37r‘^D 

N ^2'KrL 2(m^+mn + n^) 

The nanotori in Eig. 6(b) are based on (4,4) CNTs, yield¬ 
ing AEc/N = 0.251 eV, in very good agreement with the 
fitted constant a. The residual energy term b represents 
the local contribution from the elbow joints, which be¬ 
comes negligibly small in the large L limit. 


Height Dependence 


Next, we study a series of polygonal nanotori re¬ 
sembling double-walled CNTs with their adjacent wall 
ends connected by a lip-lip interaction^ in the form of 
a graphitic network, shown in Eig. 6(c). With increas¬ 
ing height H of the nanotorus, the curvature energy per 
atom AEc/N will be increasingly dominated by the cen¬ 
tral tubular part. In analogy to arguments used for nan¬ 
otori with varying L, the total number of atoms N is pro¬ 
portional to H except a finite number at the double-wall 
nanotube ends. Then, the curvature energy per atom can 
be expressed by AEc{H)/N = a' + b' /H. This behav¬ 
ior, with the values a' = 0.080 eV and b' = 1.658 eV-A, 
is reproduced by the dashed line in Eig. 6(c) for a par¬ 
ticular series of nanotori depicted in the Supplemental 
Material^^. 

Similar to the case discussed in Eig. 6(b), AEc/N can 
be approximated by the sum of curvature energies of two 
nanotubes of length i7, the inner one with radius Vi and 
the outer one with radius Then, 


AEc{ri,ro 

N 


AEcjn) + AEc(ro 
’ 27r(ri + ro)H 
\/S'k^D 



Eor the nanotori presented in Eig. 6(c), the chiral indices 
are {rrii^ni) = (5,5) and {mo^rio) = (10,10), which gives 
the value AEc/N = 0.080 eV. This value, again, is in 
perfect agreement with the fitted value a' above. Again, 
the constant b' term describes the residual energy of the 
end joints connecting the inner and the outer tube. 

In the above energy estimate, we have ignored the 
inter-layer interaction between the outer and inner wall. 
If the separation Vq — Vi between the walls were as small 
as in graphite, this stabilizing interaction would reduce 
the strain energy by ?^0.03 eV/atom. In reality, a much 
smaller effect of this interaction is expected, since the 
strain at the end junctions tends to keep Vo — ri large and 
since the interaction should depend inversely on (r^—r^)^. 
In any case, the inter-wall interaction is negligibly small 
in comparison to values shown in Eig. 6(c). In other, 
non-toroidal structures including helically coiled CNTs^^, 
the significance of the inter-wall interaction has to be ac¬ 
counted for on a case-by-case basis. 


Rotational Symmetry Number Dependence 

To the best of the authors’ knowledge, all polygonal 
carbon nanotori investigated in the literature had rota¬ 
tional symmetry numbers n = 5 or n = 6. This is largely 
due to the fact that a nanotorus can be constructed in 
a cut-and-paste manner from graphene^^’^^’^^ while pre¬ 
serving its original hexagonal symmetry. Also the ma¬ 
jority of nanotori studied here has been constructed us¬ 
ing a scheme generating isomers with either Dnd oi* Dnh 
point group symmetry. This constraint reduces greatly 
the number of possible isomers, but provides for an easy 
way to characterize individual polygonal nanotori by the 
relative positions of non-hexagonal rings within a rota¬ 
tional unit cell. With only selected results for n = 5 and 
n = 6 at hand, no conclusions are possible regarding the 
dependence of the stability on the rotational symmetry 
number. 

Such results are presented in Eig. 6(d), where we dis¬ 
play the curvature energy as a function of the rota¬ 
tional symmetry number, with n changing from n = 4 
to n = 13, while all other shape parameters are fixed. 
We focus on two structural nanotori families. The first 
family, represented by black squares in Eig. 6(d), was 
studied more extensively in the literature^^’^^’^^’^^ and 
is found to be most stable for n = 6. The second fam¬ 
ily, represented by the blue dots, consists of a series of 
nanotori with a different distribution of non-hexagonal 
rings and is energetically optimal for n = 9 instead. The 
specific shapes of these two families are presented in the 
Supplemental Material^^. 

Among the nanotori investigated in this study and 
shown in the Supplemental material^^, we find that a 
significant fraction minimized the curvature energy per 
atom for a non-hexagonal rotational symmetry. As a 
matter of fact, the distribution of ‘optimum rotational 
symmetry numbers riopt/ displayed in the Supplemen- 



tal Material^^, is roughly a Gaussian centered at Uopt^T^ 
with some extreme outliers at Uopt = 12. We found that 
the deviation of riopt from the expected value riopt = 6 
is an artifact caused by the parametrized Keating force 
field, which is penalizing bond length deviations from 
1.42 A more than bond angle deviations from 120°, while 
underestimating the penalty for out-of-plane bending and 
thus somewhat distorting the optimum geometry. 

The optimality of the n = 6 rotational symmetry num¬ 
ber can be restored in the large-torus limit as follows. 
Let us consider a cut-and-paste model of a graphitic 
torus with n = 6 in the shape displayed in Fig. 6(a). 
It is clear that all non-hexagonal rings, which deter¬ 
mine the angle ip = 120° between adjacent nanotube 
segments, will be located in the elbow joints colored in 
black, whereas the straight grey-colored segments will 
contain only hexagons. In large tori, the relative role of 
the elbow joints will play an ever diminishing role, and 
strain energy will be determined by the shape of the nan¬ 
otube segments. Deviation from n = 6 will mean that the 
nanotube segments need to be bent, which causes extra 
strain. 


CONCLUSIONS 

We have presented comprehensive analysis on the elas¬ 
tic energy of carbon nanotori containing either only 
hexagons (polyhex nanotori) or also other polygons 
(polygonal nanotori) on the basis of continuum elasticity 
theory. In polyhex nanotori, we found that depending on 
the ratio between the major radius R and minor radius 
r of the torus, the in-plane and the out-of-plane contri¬ 
butions to the total elastic energy vary significantly. The 
wide CNT rings resembling nanotori, which have been 
observed experimentally, display only negligible in-plane 
strain, whereas the in-plane strain should exceed the out- 
of plane strain for i?/r<20. In the polygonal nanotori 
studied here, the in-plane strain is rather small and the 
elastic energy obtained with the continuum method is 
shown to agree quantitatively with the results of ab ini¬ 
tio DFT calculations. We also show that the analytical 
expression Eq. (11) can serve as a quick and qualitative 
reference for the elastic energy of nanotori once the shape 


parameters (i?, r) are known. 

The capability of the current methodology is further 
demonstrated by a detailed analysis of the distributions 
of local excess energy in individual nanotori. Depending 
on the relative loci of non-hexagonal rings, the distri¬ 
bution of Gaussian curvature, mean curvature, and the 
local curvature energy can be either localized or evenly 
distributed across the nanotorus surface. This analysis 
can be extended to other 2D systems with different chem¬ 
ical composition and shape, and can also be related to the 
local stability and chemical reactivity index of different 
sites. 

We have furthermore studied three different sets of 
polygonal nanotori with varying shape parameters, in¬ 
cluding the lateral and the axial dimension of the nan¬ 
otori and the rotational symmetry number. Gontrary to 
the common perception that the most stable nanotori all 
have a six-fold symmetry, we find that in smaller polyg¬ 
onal nanotori, the optimal rotational symmetry number 
covers a wide range 4<nopt<12. Only in the large-size 
limit, when the number of non-hexagonal defects is fixed, 
n = 6 emerges as the optimum rotational symmetry num¬ 
ber. Asymptotic analysis on the variation of the other 
two shape parameters agrees quantitatively with the nu¬ 
merical results. This confirms that the current method¬ 
ology, at least for the systems investigated here, is ap¬ 
plicable across a wide length scale: From small nanotori 
where ab initio calculations are available, to mesoscopic 
tori, where continuum elasticity theory applies. Owing 
to the broad applicability, we believe that our approach 
will provide valuable results pertaining to the thermody¬ 
namical behavior of other large, experimentally observed 
carbon nanostructures, where atomic-scale treatment by 
ab initio techniques is not practical. 

ACKNOWLEDGMENTS 

This work was partly funded by the National Science 
Foundation Gooperative Agreement #EEG-0832785, ti¬ 
tled “NSEG: Genter for High-rate Nanomanufacturing”. 
Gomputational resources for this project were provided 
by the Michigan State University High-Performance 
Gomputer Genter. 


* tomanek@pa.msu.edu 

^ H. W. Kioto, J. R. Heath, S. C. O’Brien, R. F. Curl, and 
R. E. Smalley, Nature 318, 162 (1985). 

^ S. lijima, Nature 354, 56 (1991). 

^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, 1. V. Grigorieva, and A. A. 
Firsov, Science 306, 666 (2004). 

^ D. Tomanek, Guide Through the Nanoearbon Jungle^ 2053- 
2571 (lOP Publishing, Bristol, UK, 2014). 

^ M. S. Dresselhaus, G. Dresselhaus, and P. Avouris, Carbon 
Nanotubes: Synthesis, Strueture, Properties, and Appliea- 


tions (Springer, Berlin, Germany, 2001). 

® L. Liu and J. Zhao, in Syntheses and Applieations of Car¬ 
bon Nanotubes and Their Composites, edited by S. Suzuki 
(InTech, Rijeka, Croatia, 2013). 

^ M. Zheng and C. Ke, Small 6, 1647 (2010). 

® L. Chen, H. Wang, J. Xu, X. Shen, L. Yao, L. Zhu, Z. Zeng, 
H. Zhang, and H. Chen, J. Am. Chem. Soc. 133, 9654 
( 2011 ). 

® L. Liu, G. Y. Guo, C. S. Jayanthi, and S. Y. Wu, Phys. 
Rev. Lett. 88, 217206 (2002). 

M.-F. Lin, J. Phys. Soc. Jpn. 67, 1094 (1998). 



9 


M. F. Lin and D. S. Chun, Phys. Rev. B 57, 6731 (1998). 

C. G. Rocha, M. Pacheco, Z. Barticevic, and A. Latge, 
Phys. Rev. B 70 , 233402 (2004). 

A. Ceulemans, L. F. Chibotaru, and P. W. Fowler, Phys. 
Rev. Lett. 80 , 1861 (1998). 

T. A. Hilder and J. M. Hill, Phys. Rev. B 75, 125415 
(2007). 

M. T. Lusk and N. Hamm, Phys. Rev. B 76 , 125422 (2007). 

N. Hamada, S.-i. Sawada, and A. Oshiyama, Phys. Rev. 
Lett. 68, 1579 (1992). 

J. Guan, Z. Jin, Z. Zhu, G. Ghuang, B.-Y. Jin, and 

D. Tomanek, Phys. Rev. B 90 , 245403 (2014). 

L. D. Landau, L. P. Pitaevskii, A. M. Kosevich, and E. M. 
Lifshitz, Theory of Elasticity, Third Edition: Volume 7 
(Course of Theoretical Physics) (Pergamon Press, London, 
England, 1959). 

P. N. Keating, Phys. Rev. 145, 637 (1966). 

E. Artacho, E. Anglada, O. Dieguez, J. D. Gale, A. Garcia, 
J. Junquera, R. M. Martin, P. Ordejon, J. M. Pruneda, 
D. Sanchez-Portal, and J. M. Soler, J. Phys. Gond. Mat. 
20 , 064208 (2008). 

D. M. Geperley and B. J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 

J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

N. Troullier and J. L. Martins, Phys. Rev. B 43 , 1993 
(1991). 

H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13 , 5188 
(1976). 

M. R. Hestenes and E. Stiefel, J. Res. Natl. Bur. Stand. 
49 , 409 (1952). 

G. Overney, W. Zhong, and D. Tomanek, Z. Phys D: 
Atoms, Molecules and Glusters 27 , 93 (1993). 

G. Q. Ru, Phys. Rev. B 62 , 9973 (2000). 

G. Ghuang, Y.-G. Pan, and B.-Y. Jin, J. Ghem. Inf. Model. 
49 , 361 (2009). 

O. Hod, E. Rabani, and R. Baer, Phys. Rev. B 67, 195408 
(2003). 

M. Ahlskog, E. Seynaeve, R. Vullers, G. Van Haesendonck, 
A. Eonseca, K. Hernadi, and J. B Nagy, Ghem. Phys. Lett. 
300 , 202 (1999). 

R. Martel, H. R. Shea, and P. Avouris, Nature 398, 299 
(1999). 


L. Song, L. J. Gi, L. E. Sun, G. Jin, L. Liu, W. Ma, D. Liu, 

X. Zhao, S. Luo, Z. Zhang, Y. Xiang, J. Zhou, W. Zhou, 

Y. Ding, Z. L. Wang, and S. Xie, Advanced Materials 18 , 
1817 (2006). 

G. R. Ghang, L. H. Lu, J. H. Liu, and W. Ghen, Ghem. 
Phys. 393 , 123 (2012). 

J. E. Avron and J. Berger, Phys. Rev. A 51 , 1146 (1995). 
J. Berger and J. E. Avron, J. Ghem. Soc., Earaday Trans. 
91 , 4037 (1995). 

R. Martel, H. R. Shea, and P. Avouris, J. Phys. Ghem. B 
103 , 7551 (1999). 

See Supplemental Material at http://link.aps.org/ 
supplemental/10.1103/PhysRevB.000.000000 for optimum 
geometries, local Gaussian curvature and local strain dis¬ 
tribution in nanotori with different shape parameters. 

B. 1. Dunlap, Phys. Rev. B 50, 8134 (1994). 

B. 1. Dunlap, Phys. Rev. B 46, 1933 (1992). 

S. Itoh, S. Ihara, and J.-I. Kitakami, Phys. Rev. B 47, 
1703 (1993). 

S. Ihara, S. Itoh, and J.-I. Kitakami, Phys. Rev. B 47, 
12908 (1993). 

S. Itoh and S. Ihara, Phys. Rev. B 48, 8323 (1993). 

S. Itoh and S. Ihara, Phys. Rev. B 49, 13970 (1994). 

S. Ihara, S. Itoh, K. Akagi, R. Tamura, and M. Tsukada, 
Phys. Rev. B 54, 14713 (1996). 

R. Setton and N. Setton, Garbon 35 , 497 (1997). 

1. Laszlo and A. Rassat, J. Ghem. Inf. Gomput. Sci. 43, 
519 (2003). 

G. Ghuang, Y.-G. Pan, and B.-Y. Jin, J. Ghem. Inf. Model. 
49 , 1679 (2009). 

G. Ghuang, Y.-G. Ean, and B.-Y. Jin, J. Mol. Struct. 
1008 , 1 ( 2012 ). 

R. Tamura, M. Ikuta, T. Hirahara, and M. Tsukada, Phys. 
Rev. B 71 , 045418 (2005). 

E. Beuerle, G. Herrmann, A. G. Whalley, G. Valente, 
A. Gamburd, M. A. Ratner, and J. E. Stoddart, Ghem- 
istry: A European Journal 17 , 3868 (2011). 

M. T. Lusk and N. Hamm, Phys. Rev. B 76 , 125422 (2007). 
A. A. Pacheco Sanjuan, M. Mehboudi, E. O. Harriss, 

H. Terrones, and S. Barraza-Lopez, AGS Nano 8, 1136 
(2014). 



